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摘 要 : 为 研究 不 同 气 磷 水 平 下 冬小麦 将 粒 蛋 白质 含量 高 光谱 遥感 监测 模型 ， 提 高 模型 精度 ， 本 文通 过 连续 5 


Hr 


量 ， 以 相关 、 回 归 等 


和 穗 期 、 灌 浆 期 、 成 熟 期 基 


x 
定位 试验 研究 不 同 气 磷 耦合 水 平 下 ,不 同 生育 时 期 冬小麦 冠 层 光谱 反射 率 、 植 株 所 含量 以 及 成 熟 期 籽粒 蛋 


统计 分 析 方法 ,建立 基于 不 同 生 育 时 期 植株 所 含量 的 将 粒 蛋 白质 含量 监测 模型 ; A 


灰色 关联 度 分 析 ， 筛选 植株 所 含量 的 最 佳 植被 指数 ， 以 偏 最 小 二 乘 回归 法 ,建立 基于 植被 指数 的 植株 
量 监测 模型 ; 最 后 以 植株 氮 含 量 为 链接 点 ,按照 “植被 指数 一 植株 氮 含 量 一 将 粒 蛋 白质 含量 "之 间 的 联系 ， 
合 植被 指数 与 植株 氮 含 量 的 冬小麦 成 熟 期 将 粒 蛋 白质 含量 监测 模型 。 结 果 表明 : 在 拔节 期 、 孕 穗 期 、 


于 植株 气 含 量 建立 的 成 熟 期 籽粒 蛋白 质 含量 监测 模型 ， 具 有 较 好 的 监测 精度 ; 拔节 


有 、 孕穗 期 、 抽 德 期 、 灌 浆 期 、 成 熟 期 分 别 基 于 修正 叶绿素 吸收 反射 率 指数 (MCARII)、 归 一 化 差 值 叶 绿 素 指 


数 (NDCIT)、 修 正 归 一 化 差异 指数 (mNDVI)、MCARI1、NDCI 植被 指数 建立 植株 氮 含 量 监测 模型 ， 监 测 精度 (R?) 
分 别 为 0.826、0.854、0.867、0.859 和 0.819; 以 植株 氮 含 量 为 链接 点 ， 通 过 “植被 指数 一 植株 氮 含 量 一 将 粒 蛋 


粒 蛋 白质 含量 监测 模型 , R? 


2.221 gkg !、 1.825 gkg '. 


白质 含量 ”的 间接 联系 ,建立 基于 拨 节 期 、 孕 穗 期 、 抽 惩 期 、 灌 浆 期 、 成 熟 期 植被 指数 且 融 合 植株 氮 含 量 的 将 


分 别 为 0.935、0.972、0.990、0.979 和 0.936; 以 独立 数据 对 模型 进行 验证 ,模型 预 


测 值 与 实测 值 间 相对 误差 (RE) 分 别 为 11.26%、10.74%、8.41%、10.25% 和 和 11.36%， 均 方 根 误差 (RMSE) 分 别 为 


1.214 gkg. 1.767 gkg ! 和 2.137 g:-kg 1!。 说 明基 于 不 同 生 育 时 期 植被 指数 链接 植 


株 氮 含 量 可 以 对 成 熟 期 籽粒 蛋白 质 仓 量 进行 有 效 监测 ， 且 模型 具有 和 较 好 的 年 度 间 重演 性 和 品种 间 送 应 性 。 
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Abstract: Hyperspectral remote sensing data have strong band continuity, high spectral resolution and rich spectrum 
information. It can rapidly and nondestructively acquire vegetation information and it is an reliable real-time technology 
applicable in the monitoring and management of crop growth. Grain protein content (GPC) is an important indicator for wheat 
quality. Early detection of GPC of wheat using hyperspectral remote sensing data can enhance effective decision-making to 
support the acquisition and processing of high quality wheat. The objectives of this study were to establish GPC estimation 
model based on winter wheat canopy hyperspectral reflectance at different growth stages with different rates of nitrogen or 
phosphorus applications. The overall goal was to improve forecast precision of GPC estimation model at different growth stages of 
winter wheat. Thus experiments were carried out in 2009-2014 at Northwest A&F University, Shaanxi Province. The treatments 
included different winter wheat varieties with various drought resistances under five nitrogen fertilizer application rates (0, 75, 150, 
225 and 300 kg:hm ? pure nitrogen) and four phosphorus application rates (0, 60, 120 and 180 kg:hm ? P5O.). Plant nitrogen content 
(PNC) and canopy hyperspectral reflectance of different wheat cultivars were measured at jointing, booting, heading, filling and ma- 
turity stages. Also GPC was measured at maturity stage. The relationship among PNC, canopy hyperspectral reflectance and GPC 
was explored using correlation analysis, regression analysis, grey relation analysis or partial least squares. The GPC monitoring 
model was built according to the relation of "vegetation index based on canopy hyperspectral reflectance (VI)— PNC— GPC" with 
PNC as the linking point. The results showed a higher GPC prediction accuracy by GPC monitoring model based on PNC at 
jointing, booting, heading, grain-filling and maturity stages. The monitoring models of PNC at jointing, booting, heading, 
filling, maturity stages of winter wheat respectively based on modified chlorophyll absorption reflectance index (MCARI;,), 
normalized difference chlorophyll index (NDCI), modified normalized difference vegetation index (mNDVI), MCARI, and 
NDCI had better estimations of PNC, with determination coefficients (R°) of 0.826, 0.854, 0.867, 0.859 and 0.819, accordingly. 
When linked with PNC, by using the “VI— PNC— GP” method, the GPC monitoring models for the maturity stage consisted of 
combinations of VI and PNC at jointing, booting, heading, filling, maturity stages had the determination coefficients (R°) of 
0.935, 0.972, 0.990, 0.979 and 0.936, respectively. Then validatation of the models with measured values was conducted to 
verify the reliability and applicability of the models. The results showed that the relative errors (RE) between the measured and 
predicted values for the five vegetation indices were 11.26%, 10.74%, 8.41%, 10.25% and 11.36%, respectively. Then the 
corresponding root mean square errors (RMSE) were 2.221 g-kg !, 1.825 gkg™', 1.214 g-kg !, 1.767 g-kg and 2.137 g-kg !. It 
therefore suggested that MCARI;, NDCI, mNDVI, MCARI' and NDCI vegetation indices were the most suitable model for 
monitoring winter wheat GPC at jointing, booting, heading, filling and maturity stages, respectively. There was higher 
prediction precision with different vegetation indices at different growth stages monitoring winter wheat GPC under different 
N and P rates. Furthermore, the monitoring model based on different vegetation indices at different growth stages had a higher 
prediction accuracy. This results provided technical support for GPC monitoring of winter wheat at different fertilization and 
different growth stages. 

Keywords: Winter wheat; Hyperspectral remote sensing; Canopy spectral reflectance; Grain protein content; Plant nitrogen 


content; Vegetation index 


籽粒 蛋白 质 含量 (grain protein content, GPC)zÉ 
Efi] NAE Triticum astivum) 品 质 的 主要 因素 之 一 。 发 
展 优质 专用 小 麦 ， 在 保障 粮食 安全 的 同时 提高 粮食 
品质 ， 是 提高 农业 生产 效益 的 必然 趋势 "1。 传 统 小 
麦 籽 粒 品质 分 析 方 法 难以 大 范围 实时 监测 小 麦 籽粒 
蛋白 质 含 量 ， 且 费时 费 工 。 因此 ， 亚 需 一 种 在 收获 前 
对 籽粒 蛋白 质 含量 大 面积 实时 无 损 监测 的 新 技术 。 
高 光谱 遥感 凭借 波段 连续 性 强 、 e 
势 , 为 实时 快速 监测 作物 营养 状况 、 群 体 长 势 、 
量 及 品质 等 领域 提供 有 效 途径 :"。 

已 有 科学 家 通过 遥感 技术 实现 了 作物 GPC 的 监 
测 研究 ， 并 证 实 其 可 行 性 。Hansen 等 提出 以 冠 层 
光谱 反射 率 和 偏 最 小 二 乘法 预测 小 麦 GPC。Basnet 
等 中 研究 表明 基于 卫星 影像 光谱 信息 监测 谷类 作物 
GPC 具有 和 良好 效果 。Reyniers 等 ("基于 归 一 化 差 值 


植被 指数 (normalized difference vegetation index, 
NDVD 实 现 了 冬小麦 GPC 的 遥感 监测 。Liu SU 
于 Landsat TM 影像 构建 结构 不 敏感 指数 (structure 
insensitive pigment index,，SIPD 在 灌浆 期 成 功 预测 小 
Æ GPC, Wang 等 中 "研究 报道 小 麦 开 花期 植物 色素 
比率 指数 (plant pigment ratio, PPR) 与 成 熟 期 GPC 有 
较 好 的 相关 性 。Zhao 等 "基于 氮 素 运转 原理 , 指出 
冬小麦 开花 期 水 分 指数 (water index，WD 与 GPC R 
有 和 较 好 的 相关 性 。Pettersson 等 0 指出 可 以 利用 转化 
叶绿素 吸收 反射 指数 (transition chlorophyll absorp- 
tion reflection index, TCARD 对 大 麦 籽粒 粗 蛋 白 含量 
进行 预测 。 

薛 利 红 等 上 研究 认为 拔节 期 至 灌浆 期 水 稻 
(Oryza sativa) 冠 层 光 谱 反 射 率 与 成 熟 期 GPC 的 相关 
系数 达 0.80 以 上 (P<0.01)。 王 纪 华 等 中 提出 通过 
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侍 革 867 


820-1 140 nm 波段 冠 层 光谱 反射 率 预测 叶片 合 氮 量 
间接 预测 GPC 的 设想 。 田 永超 等 "4 研究 指出 小 麦 
抽穗 后 比值 植被 指数 (RVIisoosio) 与 GPC 呈 极 显著 指 
数 关系 。 薛 利 红 等 上 通过 冠 层 光谱 反射 率 监测 叶片 
氮 合 量 间 接 预测 小 麦 GPC 预测 精度 达 0.88 FPE 
等 (基于 TM 遥感 影像 构建 籽粒 氮 运 转动 态 的 小 麦 
GPC 监测 模型 ， 均 方 根 误差 幅度 为 0.47%~0.59%。 金 
秀和 良 等 5 将 灰色 关联 度 分 析 (grey relational analysis, 
GRA) 与 偏 最 小 二 乘 回归 法 (partial least squares, PLS) 
相 结合 ， 通 过 比率 指数 与 乘积 指数 改进 小 麦 GPC 的 
估算 方法 。 李 振 海 等 名 基于 氮 素 运转 理论 ， 整合 
GRA-PLS， 耦 合 叶片 所 含量 监测 模型 与 籽粒 氮 素 运 
转 模 型 ， 建 立 小 麦 GPC 监测 模型 。 

综 上 所 述 ， 目 前 利用 高 光谱 遥感 技术 监测 GPC, 
多 基于 某 一 关键 生育 时 期 或 植被 指数 与 GPC 的 相关 
性 直接 监测 ,但 这 些 研究 对 不 同 生育 时 期 植株 氮 素 
状况 与 植被 指数 间 关 联 程度 的 考虑 较 少 ， 导 致 模型 
的 稳定 性 和 重复 性 较 差 。 因 此 ， 本 文 连续 5 年 定位 
研究 不 同 氮 磷 水 平 下 ,不 同 生 育 时 期 冬小麦 冠 层 光 
谱 反 射 率 、 植 株 氮 含 量 (plant nitrogen content, PNC) 
及 成 熟 期 GPC, 首先 建立 基于 不 同 生 育 时 期 植株 氮 
合 量 的 GPC 监测 模型 ; 然后 分 析 不 同 生育 时 期 植被 
HAS PNC 的 灰色 关联 度 ， 筛 选 最 佳 植被 指数 ， 通 
过 偏 最 小 二 乘 回归 方法 ， 建 立 基 于 最 佳 植被 指数 的 
不 同 生 育 时 期 PNC 监测 模型 ; 最 后 以 PNC 为 链接 点 ， 
按照 “植被 指数 -植株 氮 含 量 -籽粒 蛋白 质 合 量 ” 之 间 
的 联系 ,建立 既 包 含 植被 指数 又 融合 不 同 生 育 时 期 
PNC 的 冬小麦 成 熟 期 GPC 监测 模型 ， 并 以 独立 数据 
对 模型 进行 验证 , 增强 模型 的 农学 解释 机 理 , 为 高 


光谱 遥感 技术 监测 小 麦 品质 提供 理论 依据 和 技术 支 
撑 ， 促 进 冬 小 麦 品质 管理 的 动态 调控 。 
1 材料 与 方法 
1.1 试验 地 概况 

试验 设 于 西北 农林 科技 大 学 北 校区 ，108?*10 了 E， 
34? 10'N, 海拔 454.8 m, 温带 大 陆 性 季风 型 气候 , 年 平均 
温度 12-14'C, 无 霜 期 220 d, FRÆ 1 400 mm, FF 
均 降 水 量 621.6 mmo 试验 地 为 粉 砂 黏 塘 土 , 0-20 cm EE 
土壤 养分 含量 为 : 有 机 质 14.26 gkg', 全 和 氮 0.90 gkg 
碱 解 氮 36.00 mgkg !, 速效 磷 17.64 mg'kg'; 20-40 cm E 
层 养分 合 量 为 : 有 机 质 10.04 gkg ', 全 氮 0.62 gkg”, 
碱 解 氮 25.29 mg-kg !， 速 效 磷 23.85 mg'kg ^c 
12 ”试验 设计 

试验 于 2009 一 2014 年 实施 , 采取 随机 区 组 设计 。 选 
择 当 地 已 有 较 大 推广 面积 的 冬小麦 品种 ,依据 不 同 品种 
的 相对 抗旱 指数 ,每 年 设置 2 个 不 同 抗旱 程度 的 品种 。 
依据 该 地 区 小 麦 生产 措施 , 设置 5 个 氮 素 水 平 , 4 个 磷肥 
水 平 , 3 组 重复 。 小 区 面积 为 3 mx10 mo 5 个 氮肥 (尿素 ， 
46% NN) 水 平分 别 为 : No( 不 施肥 )、N1(75 kg:hm ? £i 
氮 )、N2(150 kg:hm” 纯 氮 )、N3(225 kg:hm” 纯 氮 ) 和 
N4(300 kg-hm” 纯 氮 )， 总 氮肥 60% 作 为 基肥 ，40% 作 
为 追肥 ， 于 返青 后 拔节 前 施 入 。4 个 磷肥 (16% 
CaP;H4O8) 水 平分 别 为 : Po( 不 施肥 )、P1(60 kghm 
P,O;. P.,(120 kg'hm” P205) P3(180 kg-hm ? P50,), 
磷肥 作为 底肥 一 次 性 施 入 。 不 施 钾肥 ， 其 他 田间 管 
理 同 黄土 高 原 高 产 田 。 播 种 量 均 为 187.50 khm”, 
每 小 区 播种 量 约 0.56 kg。 不 同年 份 供 试 品种 、 播 种 、 
收获 及 数据 采集 日 期 见 表 1。 


表 1 不 同年 份 田间 试验 基本 情况 


Table 1 Basic information of the experiment in different years 
年 份 品种 Variety 播种 日 期 (年 -月 -日 ) 收获 日 期 (年 -月 -日 ) 数据 采集 日 期 (月 -日 ) 
Nar 强 抗 旱 型 弱 抗 旱 型 Sowing date Harvesting date Data collection time 
Strong drought resistant Weak drought resistant (year-month-day) (year-month-day) (month-day) 

2009—2010 EE 47 JM47 西 农 979 XN979 2009-10-09 2010-06-05 04-07, 04-22, 05-07, 05-20, 05-27 
2010—2011 EXE 47JMA7 西 农 979 XN979 2010-10-13 2011-06-10 04-09, 04-23, 05-07, 05-18, 06-03 
2011—2012 普 冰 9946 PB9946 西 农 2000 XN2000 2011-10-16 2012-06-08 04-08, 04-23, 05-05, 05-19, 06-03 
2012—2013 周 麦 18 ZM18 小 优 22 XY22 2012-10-10 2013-06-09 04-03, 04-17, 05-03, 05-18, 05-30 
2013—2014 长 旱 58 CHs8 BEE LM9 2013-10-08 2013-06-11 04-05, 04-20, 05-04, 05-19, 06-04 


1.3 ”测定 指标 及 方法 
131 ” 冠 层 光谱 反射 率 测量 

冬小麦 冠 层 光谱 反射 率 采 集 用 美国 Analytical 
Spectral Devices (ASD) 公 司 Field-spec Pro FR-2500 型 
背 挂 式 野外 高 光谱 辐射 测量 仪 ， 波 段 为 330~2 500 nm, 
350-1 000 nm 光谱 采样 间隔 为 1.4 nm， 光 谱 分 辨 率 


为 3 nm; 1 000~2 500 nm 光谱 采样 间隔 为 2 nm， 光 
谱 分 辨 率 为 10 nm。 分别 于 冬小麦 拔节 期 、 孕 穗 期 、 
抽穗 期 、 灌 浆 期 和 成 熟 期 ， 天 和 气 晴 朗 、 无 风 时 测量 
冠 层 光谱 反射 率 ， 适宜 时 间 为 10:00 一 14:00。 测 量 
时 传感器 探头 垂直 向 下 ， 光 谱 仪 视 场 角 为 7.5°, EE 
冠 层 顶 垂 直 高 度 约 25-30 cm， 地 面 视 场 范 围 直径 
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为 0.5 m。 每 个 观测 点 记录 10 个 采样 光谱 ,每 小 区 
测 3 次 重复 ， 取 算术 平均 值 作为 该 观测 点 冠 层 光谱 
反射 率 。 测 量 过 程 中 ,于 每 一 处 理 测量 完毕 后 ， 太 
时 进行 标准 白板 校正 (标准 白板 反射 率 为 1， 因 此 所 
得 目标 物 光谱 为 相对 反射 率 )， 以 便于 准确 测量 下 
一 处 理 。 
1.3.2 ”植株 氮 含量 (PNC) 测 定 

与 冠 层 光谱 反射 率 测定 同步 进行 。 测 完 光 谱 反 
射 率 后 ,在 同一 取样 范围 内 选取 长 势 一 致 具有 代表 
性 的 植株 20 株 ， 将 不 同 植株 按 器 官 分离 , 在 105 'C 
杀青 30 min, 85 C 烘 干 至 恒 重 ,将 不 同 器 官 分 别称 
重 后 粉碎 ，H2SO4s-H20, 消 煮 ,用 瑞士 Buchi 全 自动 


B-339 凯 氏 定 氮 仪 测定 PNC。 
1.3.3 ”籽粒 蛋白 质 合 量 (GPC) 测 定 

将 成 熟 期 测 产 小 麦 风 干 脱粒 ， 粉碎 后 以 
HSO4-H2O; 消 煮 ， 用 瑞士 Buchi 全 自动 B-339 凯 氏 
定 氮 仪 测定 。 小 麦 GPC(grain protein content, g-kg !) 
是 籽粒 全 氮 含 量 与 5.823 乘积 请。 
14 植被 指数 选择 

本 文 在 综合 前 人 关于 植被 指数 、 植 株 氮 含 量 、 
籽粒 蛋白 质 含量 研究 的 基础 上 ， 对 氮 素 营养 相关 植 
被 指数 进行 归 类 筛选 , 并 应 用 于 不 同 生 育 时 期 对 籽 
粒 蛋白 质 含量 监测 建 模 中 。 本 文选 择 主要 植被 指数 
计算 方式 及 来 源 见 表 2. 


表 2 本 研究 采用 的 高 光谱 植被 指数 


Table 2  Hyperspectral vegetation indices used in this study 


植被 指数 计算 公式 来 源 
Vegetation index Formula Reference 
修正 归 一 化 差异 指数 Modified normalized difference vegetation index (mNDVI) mNDVI=(R7s0—R705)/(R7s0+ R70s—2Ra4s) 22] 
归 一 化 差 值 叶绿素 指数 Normalized difference chlorophyll index (NDCI) NDCI=(R708—R665)/(R708+R665) 23] 
归 一 化 差异 指数 Normalized difference index (NDI) NDI=(Rsoo 一 Reso)/CRsoo+R6so) 22] 
修正 三 角 植 被 指数 Modified triangle vegetation index (TVIBL) TVIBL=0.5[120(R7s0—Rs00)—200(Re70—Rss0)] 24] 
u 绿 素 —— E "m . 

修正 叶绿素 吸收 反射 率 指数 Modified chlorophyll absorption reflectance index MCARI;-71.2[2.5(Rsoo—-Re70)-1.3(Rso9—R500)] 25] 

(MCARI;) 
红 边 归 一 化 植被 指数 Red edge normalized difference vegetation index (NDVI;os) NDVLos-(R75o—R305)/ (75g R95) 22] 
归 一 化 绿 度 指 数 Green normalized difference vegetation index (GREEN-NDVI) GREEN-NDVI-(R;5o-R550)/(Rz759-Rs50) 26] 
红 边 位 置 指数 Red edge position index (REP) REP=700+40[(R670+R780)/2-R700]/(R740-R700) 27] 
修正 简单 比值 植被 指数 Modified simple ratio index (mSRI) mSRI=(R750—R445)/(R7s0tR445) 22] 
结构 不 敏感 色素 指数 Structure insensitive pigment index (SIPI) PISI=(R570—R531)/(R570+R531) 28] 
植被 色素 比率 Plant pigment ratio (PPR) PPR=(Rsso 一 R450)/(R5s0+R4s0) 29] 


及 为 光谱 反射 率 ， 其 后 数字 为 波长 。R is spectral reflectance, the number following R is wavelength. 


15 灰色 关联 度 

关联 度 是 衡量 不 同 因素 间 关 联 程度 的 度量 , 若 
因素 间 变 化 趋势 一 致 ， 则 关联 度 较 高 ; 否则 关联 度 
较 低 。 厌 色 关 联 度 是 根据 系统 内 不 同 因素 之 间 发 展 
趋势 的 相似 或 相 异 程度 ,衡量 因素 间 关 联 程度 的 系 
统 分 析 方 法 ， 是 在 反映 系统 行为 的 特征 数据 序列 和 
影响 系统 行为 的 相关 因素 序列 中 找 出 关联 性 ， 其 思 
路 是 计算 主要 系统 行为 的 特征 数据 序列 及 影响 因素 
间 的 相关 性 ， 并 对 这 些 数据 进行 处 理 ， 进 而 定量 描 
述 不 同 因素 对 系统 的 影响 "I。 如 果 将 关联 度 认 为 是 
因素 间 的 距离 ， 则 其 只 能 为 正 值 Sl。 本 文 以 PNC 与 
植被 指数 间 的 关联 度 来 评价 二 者 之 间 同 步 变 化 程 
度 。 首 先 分 析 其 相关 性 ,选择 与 PNC 呈正 相关 的 植 
被 指数 进行 灰色 关联 度 分 析 。 
1.6 ”数据 处 理 与 分 析 

将 试验 所 得 数据 ,依据 不 同年 份 、 不 同 气 磷 耦 


合 水 平 、 不 同 生 育 时 期 、 不 同 品种 将 冠 层 光谱 反射 
率 、PNC 及 GPC 分 别 汇 总 。 用 ViewSpec 软件 对 冠 
层 光 谱 反 射 率 进行 预 处 理 。 以 2012—2013 年 数据 为 
例 分 析 不 同 氮 磷 水 平 对 成 熟 期 小 麦 GPC 的 影响 ; 综 
合 2010 一 2013 年 大 田 实测 数据 ， 对 不 同 生 育 时 期 
PNC 与 成 熟 期 GPC 进行 相关 分 析 ,， 建立 基于 PNC 
的 GPC 监测 模型 ; 同时 对 不 同 生 育 时 期 冬小麦 PNC 
与 植被 指数 进行 相关 及 灰色 关联 度 分 析 ， 在 不 同 生 
育 时 期 选择 最 佳 植被 指数 , 通过 偏 最 小 二 乘 回归 方法 ， 
建立 基于 植被 指数 的 PNC 监测 模型 ， 并 以 2009—2010 
年 独立 数据 进行 精度 验证 ; 最 后 以 PNC 为 链接 点 , 按 
照 “ 植 被 指数 一 植株 氮 含 量 一 籽粒 蛋白 质 含 量 ” 的 联 
R, 将 基于 PNC 的 GPC 模型 与 基于 植被 指数 的 PNC 
模型 组 合 ， 建 立 基于 植被 指数 且 融 合 PNC 的 成 熟 期 
GPC 监 测 模型 , 并 以 2013 一 2014 年 独立 数据 对 模型 
进行 验证 。 
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白质 含量 的 高 光谱 遥感 监测 模型 六 hi sy -141869 


2 结果 与 分 析 
21 不 同 氮 磷 水 平 对 冬小麦 籽粒 蛋 
的 影响 
以 2012—2013 年 ' 周 麦 18: 和 :小 优 22 RAH GPC 
为 例 ,， 分 析 不 同 氮 磷 水 平 对 冬小麦 GPC 的 影响 (图 1)。 研 
究 表 明 : 在 相同 氮肥 水 平 下 ， 随 着 施 磷 量 增加 (N,P0、 
NP1、N2P;、N2P;3), GPC 呈 显 著 增加 趋势 (P<0.05), 如 ' 周 
麦 18'N?P;(139.22 gkg 分 别 较 N2Pu(114.33 gkg P. 
N;P,(123.81 g'kg ). N;P,(129.09 g'kg -处理 增加 
24.88 gkg . 15.40 gkg ^. 10.13 gkg, 增幅 分 别 为 
17.87%, 11.06%, 7.2896; 在 相同 磷肥 水 平 下 ， 随 着 
氮肥 供应 量 的 增加 (NoP;、NiP;、N2P;、N3P;、NaP;), 


200 
a 周 麦 18 ZM18 a 小 价 22 XY22 


白质 含量 (GPC) 


籽粒 蛋白质 含量 
Grain protein content (g-kg ') 


N2P0 N2PI N2P2 N2P3 


Fig. 1 


Nos Nis NA. Nas Nau 238079 0 kg-hm ^, 75 kg-hm ^, 150 kg-hm ^, 


GPC 亦 呈 显著 递增 趋势 (P<0.03)， 如 ' 周 麦 18" NP 
(17374 gkg ) 分 别 较  NoP(9241 gkg) 、 NP» 
(115.10 gkg ')、 N;P4(129.09 gkg ')、 N:P, (164.11 gkg ') 
增加 81.33 g'kg', 58.64 gkg |. 44.65 gkg, 9.63 g'kg', 增 
幅 为 46.8190. 33.75%, 25.70%, 5.5496. “JME 22” 
亦 呈 现 相 似 趋势 。 说 明 随 着 氮 磷 供应 量 的 增加 ， 
GPC 呈 逐 渐 递 增 趋势 ， 但 增加 幅度 逐渐 降低 ; 增 施 
氮肥 对 GPC 提高 效果 较 增 施 磷 肥 幅 度 更 加 明显 。 

其 他 年 度 间 , 不同 品种 也 表现 出 相似 趋势 。 不 同 品 
种 间 GPC 存在 一 定 差异 ,这 与 不 同 品种 对 氮 磷 肥 
的 吸收 利用 效率 有 关 ， 同时 与 作物 品种 遗传 性 有 一 
定 关系 。 


P2NO  JP2NI P2N2  P2N3 P2N4 
施肥 水 平 Fertilizer level 


1 不 同 氮 磷 水 平 下 冬小麦 籽粒 蛋白 质 (GPC) 含 量 (2012 一 2013 年 ) 


Grain protein contents at different N and P rates of wheat in 2012 and 2013 
225 kg-hm “和 300 kg:hm” 纯 氮 , Pos Pis Pas P5 23579 0 kg-hm ? 


、60 kg-hm ^, 


120 kg.hm ”和 180 kg-hm ? P2Oso 不 同 小 写字 母 表示 处 理 间 差异 显著 (P<0.05)。No, Ni, No, Ns, and N, indicate 0, 75, 150, 225 and 300 kg-hm ? pure nitrogen, 
respectively. Po, P;, P», and P; indicate 0, 60, 120 and 180 kg.hm ? P20;, respectively. Different lowercase letters mean significant differences at 0.05 probability. 


2.0 ”基于 植株 氮 合 
(GPC) 监 测 模型 
综合 2010—2013 年 大 田 实 测 数 据 (n=120)， 对 不 
同 生 育 时 期 PNC 与 成 熟 期 GPC 进行 相关 分 析 ( 表 3)。 
由 表 3 可 知 : 拔节 期 、 孕 穗 期 、 抽 穗 期 、 灌 浆 期 、 成 
熟 期 PNC 与 成 熟 期 GPC 均 达 极 显著 相关 (P<0.01), 相 
关系 数 (7) 分 别 为 0.781、0.810、0.851、0.876 和 0.849, 
表明 可 以 基于 不 同 生育 时 期 PNC 对 成 熟 期 GPC 进 
行 预测 。 经 回归 分 析 ， 建 立 基 于 不 同 生 育 时 期 PNC 


量 (CPNC) 的 籽粒 蛋白 质 含 量 


表 3 基于 不 同 生 育 时 期 冬小麦 植株 氮 含 量 (x) 的 籽粒 蛋白 质 含 


的 GPC 模型 ,不同 生育 时 期 均 以 二 项 式 范 数 (y= 
ax+bx”+c) 拟 合 效果 最 好 ， 拔 节 期 、 人 孕穗 期 、 抽穗 期 、 
灌浆 期 和 成 熟 期 拟 合 精度 (determination coefficient， 
R”) 分 别 为 0.636、0.661、0.767、0.811 和 0.770。 
拔节 期 、 孕 穗 期 拟 合 精度 较 低 ， 这 或 许 与 该 阶段 冬 
小 麦 植株 积累 养分 主要 用 于 营养 生长 ， 而 抽穗 期 
到 灌浆 期 ， 由 营养 生长 转 入 生殖 生长 ， didis 
运输 并 积累 ， 所 以 灌浆 期 PNC 5 GPC 的 拟 合 精 
更 好 。 


量 (y) 模 型 (n=120) 


Table 3 Monitoring models of grain protein content (y) of winter wheat based on plant nitrogen content (x) at different growth stages (n=120) 


生育 时 期 相关 系数 (x) 拟 合 方程 决定 系数 (R”) 
Growth stage Correlation coefficient Equation Determination coefficient 
拔节 期 Jointing 0.781” y-7-0.748x^*7.561x-2.315 0.636 
孕穗 期 Booting 0.810” y7-0.297x^--4.381x-4.658 0.661 
抽穗 期 Heading 0.851™ y7-0.922x^8.688x—3.459 0.767 
灌浆 期 Filling 0.876" y7-0.844x^*9.595x-10.014 0.811 
成 熟 期 Maturity 0.849" y7-1.0273^7.794x*2.167 0.770 


** 表 示 在 0.01 水 平 上 差异 显著 。** represents significant correlation at 0.01 level. 
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2.5 ”不同 生育 时 期 植被 指数 与 植株 氨 含 量 (PNC) 的 

灰色 关联 度 分 析 

以 2010—2013 年 大 田 实 测 数据 (n=120)， 对 不 
同 生 育 时 期 冬小麦 PNC 与 植被 指数 进行 相关 分 析 
( 表 4)。 由 表 4 可 知 ， 在 不 同 生育 时 期 除 PPR PNC 
没有 显著 相关 (P>0.05)、SIPI 与 PNC 呈 极 显著 负 相 
关 (P<0.01) 以 外 , mNDVI、NDCI、NDI、TVIBL、 
MCARI1!、NDVI0;s、GREEN-NDVI、mSRI、REP 等 
植被 指数 均 与 PNC 呈 极 显著 正 相关 (P<0.01)。 所 以 ， 
依据 关联 度 理论 ， 选 择 与 PNC 具有 正 相关 的 植被 指 
数 进行 灰色 关联 度 分 析 。 

表 5 为 不 同 生 育 时 期 植被 指数 与 PNC 的 灰色 关联 


度 分 析 (n=120)。 由 表 5 可知， 植被 指数 与 PNC IREX 
联 度 大 小 顺序 拔节 期 依次 为 MCARIL> 
NDCI»mNDVI-NDP-TVIBL»NDVL52GREEN-NDVI- 
REP-mSRI # f& RH fk XR X NDCIPMCARI;- 
mNDVI-NDI-TVIBL»2NDVL*REP-GREEN-NDVI- 
mSRI， 抽 穗 期 依次 为 mNDVI-NDCI-NDI- 
MCARI!>TVIBL>NDV1;0s>GREEN-NDVI>REP>mSR], 
灌浆 期 依次 为 MCARI>mNDVI>REP>NDI> 
TVIBL>NDCI>NDVI ps>GREEN-NDVI>mSRI, 成 熟 
期 依次 为 NDCI>mNDVI>NDI>TVIBL>MCARII> 
NDV1l0s>REP>GREEN-NDVI>mSRI。 分 别 于 不 同 生 育 
时 期 选择 与 PNC 灰色 关联 度 最 高 的 植被 指数 用 于 对 


PNC 进行 监测 。 


表 4 不 同 生 育 时 期 植被 指数 与 植株 氮 含 量 (PNC) 的 相关 性 (z=120) 
Table 4 Correlation between vegetation indexes and plant nitrogen content (PNC) of winter wheat at different 
growth stages (1-120) 


植被 指数 拔节 期 孕穗 期 抽穗 期 灌浆 期 成 熟 期 

Vegetation index Jointing Booting Heading Filling Maturity 

修正 归 一 化 差异 指数 Modified normalized difference vegetation index (mNDVI) 0.891” 0.875" 0.9317 0.884" 0.892" 

归 一 化 差 值 叶绿素 指数 Normalized difference chlorophyll index (NDCI) 0.897” 0.924” 0.904™ 0.831" 0.905" 

归 一 化 差异 指数 Normalized difference index (NDI) 0.868" 0.842" 0.892" 0.851" 0.856" 

修正 三 角 植 被 指数 Modified triangle vegetation index (TVIBL) 0.862” 0.835" 0.841" 0.847" 0.854" 

修正 叶绿素 吸收 反射 率 指数 Modified chlorophyll absorption reflectance index 0.909" 0.881" 0.889" 0.927" 0.841 
(MCARI;) 

红 边 归 一 化 植被 指数 Red edge normalized difference vegetation index (NDVl;os) 0.8537 0.833" . 0.834" 0.821" 0.781" 

归 一 化 绿 度 指 数 Green normalized difference vegetation index (GREEN-NDVI) 0.844" 0.797” ”0.829” 0.782" 0.748" 

结构 不 敏感 色素 指数 Structure insensitive pigment index (SIPI) -0.681" | -0.673"  -0.554" -0.467” | -0.616" 

修正 简单 比值 植被 指数 Modified simple ratio index (mSRI) 0.815" 0.7118" | 0.793" 0.627" 0.714" 

红 边 位 置 指数 Red edge position index (REP) 0.826" 0.832" . 0.801" 0.879" 0.751" 

植被 色素 比率 Plant pigment ratio (PPR) —0.383N5  —0.4045$  —0.363™5 —0.28555 — —0.154N$ 


* 和 和 ** 分 别 表示 在 59081 1% 水 平 上 差异 显著 , NS 表示 相关 性 不 显著 。* and ** represent significant correlation at 5% and 1% levels, respec- 


tively. NS represents no significant correlation. 


ze 5 不 同 生育 时 期 植被 指数 与 植株 氮 含 量 (PNC) 的 灰色 关联 及 大 小 顺序 (n=120) 
Table 5 Grey correlation degrees and rank of vegetation indexes for plant nitrogen content (PNC) of winter wheat at different 
growth stages (n=120) 


Jointing Booting Heading Filling Maturity 
EHE n. " m u " 
植被 指数 灰色 关联 度 灰色 关联 度 色 关联 度 灰色 关联 度 灰色 关联 度 
Vegetation index (£-0.5) 顺序 (5=0.5) 顺序 (£-0.5) 顺序 (&=0.5) 顺序 (£-0.5) 顺序 
Grey Order Grey Order Grey Order Grey Order Grey Order 
correlation correlation correlation correlation correlation 
mNDVI 0.915 3 0.898 3 0.957 1 0.908 2 0.916 2 
NDCI 0.943 2 0.950 1 0.929 2 0.852 6 0.940 1 
NDI 0.891 4 0.863 4 0.916 3 0.873 4 0.878 3 
TVIBL 0.885 5 0.856 5 0.862 5 0.869 5 0.876 4 
MCARI; 0.957 1 0.905 2 0.913 4 0.926 1 0.862 5 
NDVI7o5 0.875 6 0.854 6 0.855 6 0.841 7 0.859 6 
GREEN-NDVI 0.866 7 0.816 8 0.850 7 0.823 8 0.834 8 
mSRI 0.835 9 0.730 9 0.827 9 0.797 9 0.791 9 
REP 0.847 8 0.843 7 0.836 8 0.903 3 0.851 7 


植被 指数 栏 各 代码 的 意义 见 表 4。 Meanings of the codes of vegetation indexes are shown in the Table 4. 
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2.4 不同 生育 时 期 PNC 高 光谱 监测 模型 的 建立 与 

验证 

通过 表 5 灰色 关联 度 分 析 ， 获 得 不 同 生育 时 期 
5 PNC 稳定 相关 的 最 佳 植被 指数 ， 基 于 表 3 不 同 生 
育 时 期 PNC 与 GPC 的 相关 性 , 通过 偏 最 小 二 乘 回 归 
方法 ， 建 立 基于 最 佳 植被 指数 的 不 同 生 育 时 期 PNC 
监测 模型 (n=120)， 并 以 2009—2010 年 独立 数据 对 
模型 精度 进行 验证 (n=40)( 表 6)。 

Ex 6 拟 合 模 型 可 知 在 拔节 期 、 孕 穗 期 、 抽 穗 
期 、 灌 浆 期 、 成 熟 期 基于 MCARII、NDCI、mNDYVI、 
MCARI!、NDCI 等 植被 指数 建立 不 同 生育 时 期 冬 小 
Z PNC 模型 ， 拟 合 精度 (R3 分 别 为 0.826、0.854、 


0.867. 0.859 和 0.819， 标 准 误差 (SE) 分 别 为 0.213、 

0.191、0.136、0.177 和 0.243。 以 独立 数据 对 模型 
预测 值 进行 验证 ， 结 果 显 示 拔 节 期 、 孕 穗 期 、 抽 穗 
期 、 灌 浆 期 、 成熟 期 基于 MCARI,. NDCI, mNDVI, 
MCARII、NDCI 建立 PNC 模型 预测 值 与 实测 值 间 相 
对 误差 (RE) 分 别 为 14.08%、13.63%、10.31%、12.17% 
和 15.16%, 均 方 根 误差 (RMSE) 分 别 为 0.324 gkg, 

0.317 gkg, 0.125 g-kg ^, 0.230 gkg- 和 0.421 gkg 
表明 不 同 生育 时 期 基于 不 同 植被 指数 建立 PNC 监 
测 模型 县 有 较 好 的 拟 合 和 验证 效果 ， 且 在 不 同 生 
育 时 期 均 以 直线 方程 PNC=d.VI+e 形 式 的 拟 合 效果 


最 佳 。 


X6 不 同 生育 期 冬小麦 植株 氮 含 量 0) 与 植被 指数 (x) 的 拟 合 (n=120) 及 验证 (n=40) 
Table6 Simulating (n=120) and performance (n=40) of monitoring models between plant nitrogen content (y) and the best vegeta- 
tion index (x) of winter wheat at different growth stages 


拟 合 模型 Fitting model 


验证 模型 Performance model 


生育 时 期 植被 指数 - - 一 一 - 加 

Growth Vegetation 拟 合 方程 决定 系数 (R?) 标准 误差 (SE) — 相对 误差 (RE)  ” 均 方 根 误差 (RMSE) 
stage index Eauation Determination Standard Relative Root mean square error 

q coefficient error error (%) (gkg !) 

拔节 期 Jointing MCARII y=2.215x+4.233 0.826 0.213 14.08 0.324 

孕穗 期 Booting NDCI y=1.471x+6.483 0.854 0.191 13.63 0.317 

抽穗 期 Heading mNDVI y70.738x49.145 0.867 0.136 10.31 0.125 

灌浆 期 Filling MCARI, y-0.927x45.314 0.859 0.177 12.17 0.230 

成 熟 期 Maturity NDCI y70.898x--6.164 0.819 0.243 15.16 0.421 


植被 指数 栏 各 代码 的 意义 见 表 4。 Meanings of the codes of vegetation indexes are shown in the Table 4. 


2.5 ”籽粒 蛋白 质 含量 (GPC) 高 光谱 遥感 监测 模型 的 
建立 与 验证 

由 于 不 同 生育 时 期 PNC 与 成 熟 期 GPC 间 呈 极 
显著 相关 (P<0.01), 基于 PNC 可 以 建立 成 熟 期 GPC 
的 预测 模型 ， 基 于 不 同 生 育 时 期 PNC 拟 合 GPC 以 
y-axtbx^c 形式 效果 最 好 ; 而 基于 植被 指数 建立 不 
同 生育 时 期 PNC 遥感 模型 ， 具 有 较 高 的 拟 合 精度 和 验 
证 效果 , 不同 生 育 时 期 PNC=dxVIte 形式 拟 合 效果 最 
佳 。 因 此 ， 可 以 通过 不 同 生育 时 期 PNC 将 植被 指数 与 
GPC 进行 链接 。 按 照 “VI_PNC 一 GPC”* 的 联系 ,以 
PNC 为 链接 点 ,将 植被 指数 和 成 熟 期 GPC 链接 ,将 基 
于 植被 指数 的 PNC 监测 模型 与 基于 PNC 的 籽粒 蛋白 
质 监测 模型 组 合 ,建立 基于 植被 指数 且 融 合 不 同 生 育 
时 期 PNC 动态 变化 的 GPC 遥感 监测 模型 : 

GPC=ax(dxVI+e)+bx(dxVI+e) +c (1) 
式 中 : GPC 为 籽粒 蛋白 质 含量 , VI 为 植被 指数 , a、b、 
c、d、e 均 为 常数 。 

为 了 对 模型 的 稳定 性 和 适用 性 进行 评价 ， 以 
2013—2014 年 不 同 品种 实测 数据 对 上 述 优选 模型 进行 
验证 (n=40), 不 同 生 育 时 期 GPC 模型 预测 结果 见 表 7。 
由 表 7 可 知 : 拔节 期 、 孕 穗 期 、 抽 穗 期 、 灌 浆 期 和 成 熟 


期 实测 值 分 别 对 MCARI—PNC—GPC., NDCI— 
PNC—GPC, mNDVI—PNC—GPC, MCARI;— PNC— 
GPC, NDCI—PNC—GPC 路 线 拟 合 的 GPC 具有 较 好 的 
验证 效果 , 不 同 生育 时 期 模型 最 佳 监测 模型 分 别 为 : 
GPC=0.451x(0.572x MCARI1+0.244)+0.339x 


(0.572xMCARII+0.244)2+0.473 (R2-0.935) (2) 
GPC=0.627x(0.519xNDCI+0.327)+0.553x 


(0.519xNDCI+0.327)*+0.526 (R?-0.972) (3) 
GPC=0.486x(0.725xmNDVI+0.424)+0.611x 
(0.725xmNDVI-0.424) 0.418 (R?—0.990) (4) 
GPC-0.577x(0.691MCARI,;4*0.517)*0.838x 
(0.691XMCARI,;40.517)40.733 (R?-0.979) (5) 
GPC-0.591x(0.486X«NDCI--0.828)40.656x(0.486x 
NDCI-«0.828)4-0.917 (R?—-0.936) (6) 


以 独立 数据 对 模型 进行 验证 (n=40), 不 同 生 育 时 
期 监测 模型 的 预测 值 与 田间 实测 值 间 相对 误差 (RE) 
分 别 为 11.26%、10.74%、8.41%、10.25%、11.36%, 均 
方 根 误差 (RMSE) 分 别 为 2.221 g-kg . 1.825 gkg” 
1.214 gkg . 1.767 gkg ^. 2.137 gkg :。 表 明 可 以 
利用 不 同 生育 时 期 关键 植被 指数 对 不 同 氮 磷 营 养 水 
平 下 不 同类 型 冬小麦 品种 成 熟 期 籽粒 蛋白 质 进行 预 
测 ， 且 在 抽穗 期 和 灌浆 期 监测 效果 较 好 。 图 2 为 不 同 
生育 时 期 GPC 田 间 实 测 值 与 模型 预测 值 比较 结果 。 
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表 7 不 同 生育 时 期 冬小麦 籽粒 蛋白 质 


含量 监测 模型 验证 (2013 一 2014 年 ) 


Table 7 Performance of monitoring model of grain protein content of winter wheat at different growth stages (2013-2014) 


生育 时 期 技术 路 线 HEFER) 相对 误差 (RE) 均 方 根 误差 (RMSE) 
Growth stage Method Fitting precision Relative error (%) Root mean square error (g-kg !) 
拔节 期 Jointing MCARI,—PNC—GPC 0.935 11.26 2.221 
孕穗 期 Booting NDCI—PNC—GPC 0.972 10.74 1.825 
抽穗 期 Heading mNDVI—PNC—GPC 0.990 8.41 1.214 
灌浆 期 Filling MCARI,—PNC—GPC 0.979 10.25 1.767 
成 熟 期 Maturity NDCI—PNC—GPC 0.936 11.36 2.137 


PNC 为 植株 氮 含 量 , GPC 为 籽粒 蛋白 质 含量 , “技术 路 线 ” 列 的 其 他 符号 见 表 4, PNC: plant nitrogen content; GPC: grain protein content; other 


codes in the method column are shown in the Table 4. 


3 Wit 

GPC 是 评价 小 麦 品质 的 重要 指标 ， 其 气 素 积累 
主要 包括 开花 前 植株 储藏 氮 素 运转 和 花 后 氮 素 同 化 
P5882 7. BD GPC 与 不 同 生 育 时 期 植株 氮 素 营养 有 
密切 关系 ,因此 可 以 基于 不 同 生 育 时 期 PNC 对 成 熟 
期 GPC 进行 监测 所 ?9。 本 文通 过 连续 5 年 定位 研究 
不 同 生育 时 期 冬小麦 冠 层 光谱 反射 率 、PNC 以 及 成 
熟 期 GPC, 通过 相关 分 析 、 回 归 分 析 ， 建立 基于 PNC 
的 成 熟 期 GPC 监测 模型 ; 进一步 分 析 植 被 指数 与 
PNC 的 相关 性 与 灰色 关联 度 ， 筛 选 最 佳 植被 指数 ， 
以 偏 最 小 二 乘 回归 分 析 ， 建立 基 于 植被 指数 的 PNC 
监测 模型 ; 最 终 以 PNC 为 链接 点 , 将 PNC 模型 与 成 
熟 期 GPC 模型 耦合 ， 建 立 基 于 植被 指数 且 融 合 不 同 
生育 时 期 PNC 的 成 熟 期 GPC 监测 模型 ， 为 高 光谱 遥 
感 技术 监测 冬小麦 成 熟 期 GPC 增加 农学 解释 机 理 ， 
更 合理 地 阐述 了 不 同 生 育 时 期 PNC 对 成 熟 期 GPC 
形成 的 动态 理论 ,为 冬小麦 成 熟 期 GPC 快速 无 损 预 
测 提供 科学 依据 。 
3.1 基于 PNC 的 GPC 监测 模型 

不 同 生育 时 期 PNC 与 成 熟 期 GPC 均 呈 极 显著 相 
X(P«0.01), 表明 可 以 基于 不 同 生 育 时 期 PNC 实现 对 
成 熟 期 GPC 的 预测 。 这 与 前 人 研究 结果 一 致 0%?1*71。 
虽然 从 拔节 期 至 成 熟 期 均 可 依据 PNC 对 成 熟 期 
GPC 进行 预测 ,但 以 灌浆 后 期 监测 精度 较 高 。 这 是 
由 于 冬小麦 籽粒 蛋白 质 形 成 和 积累 是 动态 的 ， 在 不 
同 生育 时 期 生长 重心 不 同 ,灌浆 前 以 营养 生长 为 主 ， 
植株 氮 素 同化 积累 较 慢 ， 而 灌浆 后 , 由 营养 生长 转 
入 生殖 生长 ， 氮 素 营 养 快速 同化 并 积累 储存 于 籽粒 ， 
所 以 ,不同 生育 时 期 监测 精度 存在 差异 。 
3.2 ”基于 植被 指数 的 PNC 监测 模型 

本 文 在 对 不 同 生育 时 期 植被 指数 与 PNC 相关 
性 分 析 的 基础 上 , 分 析 植 被 指数 与 PNC 的 灰色 关联 
度 并 排序 ， 筛 选 最 佳 植被 指数 ， 这 些 植 被 指数 基于 
与 植株 氮 素 密切 相关 的 敏感 波段 构建 ， 通 过 多 波段 


组 合 增强 对 PNC 的 相关 性 ， 更 好 地 反映 不 同 生育 时 
期 PNC 状况 ; 继而 通过 偏 最 小 二 乘 回 归 建 立 基于 植 
被 指数 的 不 同 生育 时 期 PNC 监测 模型 ， 在 拔节 期 、 
孕穗 期 、 抽 穗 期 、 灌 浆 期 、 成 熟 期 基于 MCARII、 
NDCI、mNDVI、MCARII、NDCI 等 植被 指数 监测 
PNC 具有 较 好 的 监测 精度 (R*) 和 估计 标准 误差 (SE); 
同时 对 不 同 生 育 时 期 最 佳 监测 模型 以 不 同年 度 间 不 
同 品种 类 型 独立 数据 进行 评价 , 均 具 有 较 好 的 验证 
效果 ,表明 该 模型 具有 一 定 的 稳定 性 和 可 靠 性 。 
3.3 GPC 的 高 光谱 遥感 监测 模型 

冬小麦 籽粒 蛋白 质 形成 与 积累 是 一 个 动态 过 程 5? ?1， 
本 研究 以 不 同 生 育 时 期 PNC 为 切入 点 , 将 不 同 生 育 
时 期 VI 与 成 熟 期 GPC 链接 ,依据 “VI 一 PNC 一 GPC” 
的 联系 , 将 VI PNC, GPC 三 者 耦合 , 在 GPC 监测 
模型 中 引入 不 同 生 育 时 期 植株 氮 素 状况 与 VI， 并 以 
独立 数据 对 模型 稳定 性 和 适用 性 进行 评价 , 模型 具 
有 较 高 的 验证 效果 ,可 以 对 该 地 区 内 不 同 肥力 水 平 
下 不 同 抗旱 类 型 冬小麦 品种 GPC 进行 预测 。 该 方法 
以 PNC 这 一 与 GPC 密切 相关 的 农学 参数 为 中 间 变 
量 , 基于 植被 指数 建立 GPC 监测 模型 ,增强 了 GPC 
监测 模型 的 农学 解释 机 理 ， 实 现在 不 同 生 育 时 期 借 
助 高 光谱 遥感 技术 对 冬小麦 成 熟 期 GPC 的 监测 ， 同 
时 以 独立 数据 对 模型 精度 及 可 靠 性 进行 验证 , 均 有 
较 好 效果 。 
3.4 GPC 监测 模型 研究 展望 

本 文通 过 连续 5 年 对 不 同 氮 磷 水 平 下 不 同 抗旱 
类 型 的 冬小麦 品种 不 同 生 育 时 期 PNC、 冠 层 光 谱 反 


射 率 、 成 熟 期 GPC 进行 分 析 ， 基 于 “VI 一 
PNC 一 GPC” 建 立成 熟 期 GPC 监测 模型 ,具有 较 高 


的 验证 精度 。 在 建 模 过 程 中 ， 融 入 冠 层 光谱 特征 及 
植株 氮 素 营养 状况 ， 提 高 模型 的 预测 精度 和 验证 精 
度 。 由 于 受 试验 条 件 限制 ， 本 模型 的 建立 与 验证 均 
在 同一 试验 小 区 内 , 土壤 、 水 分 、 光 照 等 因素 基本 
一 致 。 今 后 模型 应 在 不 同 生 态 区 域 范围 内 ， 扩 


http://www.ecoagri.ac.cn 


第 6 期 a E$. 不 同 生 育 时 期 冬小麦 籽粒 蛋白 质 合 量 的 高 光谱 遥感 监测 模型 S c HBTS] 
n kon an : 大 建 模 因素 的 空间 差异 性 , 增强 对 模型 外 推 性 及 稳 
| K-0933 g 定性 。 另外， 本 文 建立 GPC 监测 模型 ， 仅 从 PNC 这 
Pi Pe 一 参数 出 发 ,缺乏 对 该 区 域内 常年 积温 、 降 雨量 、 
120} 土壤 状况 等 田间 生态 因素 的 考虑 。 因 此 ， 有 必要 在 
研究 生理 参数 的 同时 , 结合 区 域内 生态 因素 建立 明 
AR 确 物理 意义 的 监测 模型 ， 促 进 高 光谱 遥感 技术 在 精 
za 准 农业 中 的 应 用 。 

60 

孕 德 期 Booting 4 结论 

160|  NDCI—PNC—GPC z 

mi roms iz 1) 在 相同 氮肥 水 平 下 ，GPC 随 施 磷 量 增加 而 显 

130| RMSE=1.825 a 著 增加 (P<0.05)7.28%~17.87%; 在 相同 磷肥 水 平 下 ， 

| ZH j GPC 随 施 气量 增加 显著 递 境 5.54%~46.81% 

100} 6 (P<0.05)。 RAHA GPC 与 不 同 生 育 时 期 PNC 极 显著 

| 相关 (P<0.00)， 相 关系 数 ( 站 为 0.781~0.876， 可 以 基 

AL 于 不 同 生 育 时 期 PNC 建立 成 熟 期 GPC 模型 ， 均 以 二 

ea 项 式 画 数 0=axr+bxz+o) 拟 合 效果 最 好 ,拔节 期 、 孕 生 
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2 不 同 生育 时 期 基于 植被 指数 的 冬小麦 籽粒 蛋白 质 
含量 预测 值 与 实测 值 比较 
Fig.2 Comparison of estimated and measured grain protein 
contents in winter wheat with best vegetation index at different 
growth stages 


0.636、0.661、0.767、0.811 和 0.770。 

2) 通 过 相关 分 析 、 灰色 关联 度 分 析 ， 在 不 同 生育 
时 期 选择 最 佳 植 被 指数 ,通过 偏 最 小 二 乘 回归 方法 
建立 基于 MCARII、NDCI、mNDVI、 MCARI'、 NDCI 
的 拔节 期 、 孕 穗 期 、 抽 穗 期 、 灌 浆 期 、 成 熟 期 冬 小 
麦 PNC 高 光谱 遥感 监测 模型 ， 不 同 生 育 时 期 模型 拟 
合 精度 (R) 分 别 为 0.826、0.854、0.867、0.859 和 0.819, 
标准 误差 (SE) 分 别 为 0.213. 0.191. 0.136. 0.177 和 
0.243。 以 独立 数据 对 模型 进行 检验 ， 相对 误差 (RE) 分 
别 为 14.08%、13.63%、10.31%、12.17% 和 15.1696, 均 
方 根 误 差 (RMSE) 分 别 为 0.324 gkg’, 0.317 gkg 
0.125 gkg !、0.230 g-kg ! 410.421 gkg !。 表 明 不 同 生 
育 时 期 基于 不 同 植被 指数 建立 PNC 监测 模型 具有 较 
好 的 监测 和 验证 效果 。 

3) 以 PNC 为 链接 点 ,依据 “VI 一 PNC 一 GPC”* 的 
路 线 ， 将 基于 植被 指数 的 PNC 监测 模型 与 基于 
PNC 的 成 熟 期 GPC 监测 模型 进行 组 合 ， 建 立 基 于 
VI 且 融合 PNC 的 成 熟 期 GPC 监测 模型 ,增强 模型 
的 农学 机 理解 释 。 以 独立 数据 对 拔节 期 、 孕 穗 期 、 
抽穗 期 、 灌 浆 期 和 成 熟 期 GPC 监测 模型 进行 验证 ， 
不 同 生育 时 期 模型 预测 值 与 田间 实测 值 间 相对 误 
差 (RE) 分 别 为 : 11.269€. 10.7496. 8.41%, 10.25%, 
11.36%, 均 方 根 误差 分 别 (RMSE) 为 2.221 gkg, 
1.825 gkg”, 1214 gkg", 1.767 gkg 和 2.137 gkg™ 说 
明 可 以 通过 不 同 生育 时 期 植被 指数 监测 PNC 实现 


对 成 熟 期 GPC 的 监测 ， 且 模型 具有 和 较 好 的 稳定 性 
和 适用 性 。 
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